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SUMMARY 


This  report  presents  aerodynamic  data  on  an  NACA  0018  airfoil 
oscillating  in  pitch  in  the  NAE  2m  X  3m  Low  Speed  Tunnel.  The  purpose  of 
the  test  was  to  simulate  the  dynamic  aerodynamic  behaviour  of  a  Vertical 
Axis  Wind  Turbine  blade  section  at  the  equatorial  plane  under  the  dynamic 
stall  conditions  which  may  occur  at  low  ratios  of  tip  speed  to  wind  speed. 

Measurements  were  made  of  wing  surface  pressure  time  histories 
at  various  chordwise  locations,  during  a  complete  cycle  of  oscillation.  The 
dynamic  effects  on  nose  suction  pressures  were  seen  to  increase  their  mag¬ 
nitude  and  to  delay  flow  breakdown  to  a  higher  angle  of  attack. 

Integration  of  wing  pressures  furnished  normal  and  chordwise 
forces  during  the  angle  of  attack  cycle.  The  behaviour  of  these  parameters 
showed  that  dynamic  stall  for  wing  angles  of  attack  up  to  30  degrees, 
occurred  about  5  degrees  later  than  for  the  equivalent  steady  flow  case. 
This  phenomenon  resulted  in  an  increase  in  normal  force  of  about  20%  and 
an  increase  in  chordwise  force  of  about  40%  when  the  wing  is  pitching  to 
30  degrees. 


RESUME 

Le  present  rapport  contient  des  donnees  a6rodynamiques  sur  un 
profil  NACA  0018  oscillant  en  tangage  dans  la  soufflerie  a  basse  vitesse  de 
2  m  X  3  m  de  l’E AN.  L’essai  avait  pour  objet  de  simuler  le  comportement 
aerodynamique  dynamique  d’un  profil  de  pale  d’eolienne  a  axe  vertical 
dans  le  plan  equatorial  dans  des  conditions  de  dScrochage  dynamique,  lequel 
peut  se  produire  &  de  faibles  rapports  de  vitesse  d’extr6mit6  de  pale  sur  la 
vitesse  du  vent. 

Des  mesures  de  pression-temps  ont  ete  prises  en  differents  points 
de  la  corde  du  profil  au  cours  d’une  oscillation  complete.  Les  effets  dyna- 
miques  des  pressions  d’aspiration  s’exergant  sur  le  nez  du  profil  augmen- 
taient  alors  pour  retarder  la  deterioration  de  l’ecoulement  laminaire  aux 
angles  d’attaque  superieurs. 

L ’integration  des  pressions  s’exer^ant  sur  le  profil  a  donne  des 
forces  normales  et  qui  s’exer^aient  dans  le  sens  de  la  corde  en  fonction  de 
Tangle  d’attaque.  Ces  parametres  ont  indiqu6  que  le  d£crochage  dynamique 
pour  des  angles  d’attaque  du  profil  allant  jusqu’a  30  degr6s  se  produisait 
environ  5  degres  plus  tard  que  dans  le  cas  d’un  profil  equivalent  dont  l’ecou- 
lement  n’etait  pas  perturb^.  Ce  phenomene  s’est  traduit  par  une  augmenta¬ 
tion  d’environ  20  %  de  la  force  normale  et  de  40  %  environ  de  la  force 
s’exer^ant  dans  le  sens  de  la  corde  lorsque  le  profil  oscille  jusqu’a  30  degres 
en  tangage. 
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WIND  TUNNEL  INVESTIGATION  OF  DYNAMIC  STALL  OF  AN  NACA  0018 
AIRFOIL  OSCILLATING  IN  PITCH 


1.0  INTRODUCTION 

In  the  sphere  of  renewable  energy,  wind  power  is  considered  to  be  one  of  the  more 
appropriate  ways  of  meeting  some  of  the  future  energy  needs  of  Canada’s  remote  areas,  and  also 
those  of  underdeveloped  countries.  A  major  contribution  has  been  made  in  this  field  by  the  Low 
Speed  Aerodynamics  Laboratory  of  the  NAE,  in  the  re-discovery  and  development  of  the  Darreius 
Vertical  Axis  Wind  Turbine  (VAWT),  and  has  resulted  in  the  laboratory  playing  an  important  role 
in  the  National  Wind  Energy  Program. 

The  efforts  of  many  years  work  have  resulted  in  contributions  to  the  basic  science  of  the 
VAWT,  with  particular  reference  to  aerodynamic  performance.  Several  wind  turbine  test  installations 
are  currently  running  successfully  in  remote  parts  of  Canada,  the  largest  being  the  230  kw  and  500  kw 
machines  on  the  Magdalen  and  Prince  Edward  Islands.  The  results  of  these  experiments,  in  which 
the  benefits  of  feeding  wind-generated  electricity  into  existing  hydro  networks  as  well  as  the  assess¬ 
ment  of  the  aerodynamic  characteristics  of  a  large  machine,  have  provided  basic  research  and  design 
data  for  future  turbines.  The  development  of  an  even  larger  turbine,  on  the  south  shore  of  the 
St.  Lawrence  River  in  Quebec  is  now  underway,  and  it  is  expected  to  be  in  operation  in  1986. 

The  Magdalen  Islands  wind  turbine  is  of  the  constant  speed  type.  The  rotor,  consisting  of  two 
blades,  is  80  feet  (24.4  m)  in  diameter  at  the  equator,  and  120  feet  (36.6  m)  high.  Its  rotational  speed 
of  36  rpm  approximately,  results  in  an  equatorial  tip  speed  of  about  150  ft/sec  (45.7  m/sec).  Figure  1 
shows  a  view  of  the  turbine  and  a  typical  power  characteristic  curve  in  which  power  coefficient  Cp  is 
plotted  against  speed  ratio  Vt/V. 

The  various  operating  ranges  of  the  constant  speed  turbine  are  indicated  by  hatched  areas 
on  the  power  curve;  best  efficiency  occurs  in  the  speed  ratio  range  between  4  and  7.  In  this  operating 
range,  a  simple  momentum  theory  has  been  found  to  give  satisfactory  agreement  with  experimental 
data,  and  various  codes,  based  on  the  multiple  stream  tube  theories  of  References  (3)  and  (4)  are  now 
routinely  used  to  compute  rotor  performance.  At  lower  values  of  the  speed  ratio  (1<  Vt/V<4),  the 
turbine  is  operating  in  an  environment  where  the  angle  of  attack  of  the  blade  sections  are  at  or  above 
static  stall,  and  where  dynamic  effects,  resulting  from  cyclic  changes  of  airflow  are  at  a  maximum. 
Dynamic  aerodynamic  effects,  particularly  as  they  affect  the  airfoil  chord  force  coefficients,  are 
considered  to  have  a  significant  effect  on  turbine  maximum  power  output  and  the  extent  to  which 
this  might  be  true  is  illustrated  in  Figure  1,  in  which  dynamic  stall  effects  have  been  included  in  the 
theoretical  predictions.  At  speed  ratios  greater  than  6  the  turbine  power  diminishes  rapidly,  the  flow 
s+ate  is  considered  to  be  in  a  state  of  breakdown,  and  the  simple  axial  momentum  theories  no  longer 
apply. 


As  a  result  of  this  perceived  importance  of  dynamic  stall  effects  on  the  rotor,  an  experimental 
wind  tunnel  program  was  initiated  in  order  to  understand  more  fully  a  very  complex  phenomenon, 
and  also  to  obtain  dynamic  normal  and  chord  force  data  suitable  for  use  in  more  advanced  computa¬ 
tional  codes. 

2.0  CHARACTERISTICS  OF  UNSTEADY  FLOW  THROUGH  THE  ROTOR 

Due  to  the  cyclic  nature  of  the  relative  airflow  acting  at  a  blade  section,  fluctuations  in  both 
angle  of  attack  and  dynamic  pressure  will  occur,  resulting  in  large  variations  in  blade  sectional  load. 
The  probability  of  encountering  these  dynamic  effects  is  high  over  the  entire  range  of  operation  of 
the  rotor,  but  maximum  angles  of  attack  occur  at  all  speed  ratios  near  2  to  6  over  the  downwind 
portion  of  the  blade  cycle.  The  variation  of  angle  of  attack  over  one  revolution  and  for  various  speed 
ratios  are  shown  in  Figure  2.  At  a  speed  ratio  of  2,  the  maximum  angle  of  attack  is  30  degrees,  and 
velocity  fluctuates  by  60%.  At  a  speed  ratio  of  5,  the  maximum  angle  of  attack  is  12  degrees,  but  the 
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velocity  variation  reaches  100%.  These  calculations  have  not  included  any  secondary  effects  such 
as  those  induced  by  the  turbine  wake.  The  configuration  of  the  trailing  vortex  system  downwind 
of  the  rotor  is  cycloidal  in  form  and  its  resulting  intersection  with  the  rotating  blades  complicate  the 
aerodynamic  model,  and  may  result  in  additional  dynamic  loads.  A  water  tunnel  study  of  a  single 
blade  vertical  axis  configuration  confirmed  in  a  qualitative  way,  the  cyclic  nature  of  the  flow  through 
the  rotor,  dynamic  stall  effects,  and  a  region  of  intersection  on  the  advancing  blade  side.  This  is  also 
shown  in  Figure  2  in  which  a  large  vortex-like  disturbance  has  been  shed  into  the  flow  at  the  maxi¬ 
mum  angle  of  attack  position.  Due  to  the  low  Reynolds  number  of  this  investigation,  it  is  difficult 
to  correlate  this  disturbance  with  dynamic  stall.  The  vortex  formation  does  however  occur  once  per 
revolution,  and  suggests  that  rapid  changes  of  blade  loading  are  taking  place. 


3.0  WIND  TUNNEL  SIMULATION  OF  DYNAMIC  AERODYNAMIC  EFFECTS 

The  problems  associated  with  on-site  measurement  of  fluctuating  aerodynamic  forces  or 
pressures  on  an  operational  wind  turbine  blade  were  considered  to  be  numerous  and  probably  in¬ 
surmountable  with  current  resources  and  priorities.  A  more  acceptable  alternative  in  which  the  aero¬ 
dynamic  parameters  could  be  more  accurately  controlled,  was  a  wind  tunnel  test  in  which  the 
dynamic  aerodynamic  effects  on  a  segment  of  the  turbine  blade  could  be  measured  under  quasi- 
two-dimensional  unsteady  conditions.  Whether  or  not  this  technique  accurately  represents  the  real 
flow  at  the  blade  is  open  to  question  however,  and  the  measurement  of  actual  blade  pressures  will 
ultimately  be  required. 

Accordingly  it  was  proposed  to  oscillate  a  two-dimensional  airfoil  in  the  2m  X  3m  wind 
tunnel  in  an  angle  of  attack  time  history  which  simulates  the  local  operational  conditions  at  the 
equatorial  plane.  Surface  pressures  measured  at  a  mid -span  location  would  furnish  information  related 
to  the  chordwise  loading  on  the  blade,  and  also  the  normal  and  chordwise  forces. 

Although  dynamic  stall  on  the  wind  turbine  occurs  under  conditions  in  which  both  angle 
of  attack  and  dynamic  pressure  are  changing,  (as  in  Fig.  2)  it  would  only  be  possible  to  simulate  non¬ 
steady  angle  of  attack  variations  in  the  wind  tunnel  with  the  tunnel  speed  kept  constant. 

Referring  to  Figure  2  the  azimuthal  variation  of  angle  of  attack  and  blade  velocity  (neglecting 
induced  effects)  are  given  by  the  following  relations: 


-Sin  6 

Sin  a  =  - 

[(1  +  X2)  +  2XCos0]* 


^  =  [(1+X2)  +  2X  Cos 0 ] ,/l 


(1) 

(2) 


where  X 


,  and  Vb  is  the  resultant  flow  velocity  at  the  blade  section. 


The  dimensionless  coefficient  r  is  also  defined  which,  in  effect,  measures  the  number  of 
chord  lengths  of  air  which  have  passed  over  the  blade  during  one  revolution. 


1  T 

r  =  -  /T  Vb  •  dt 

t  O 


(3) 


where  T  is  the  period  of  one  revolution. 
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follows: 


In  terms  of  the  geometric  parameters  R/c,  Vt/V,  and  Vb/V  the  time  parameter  r  is  as 


(R/c)  a  Vb 

r=  -r-  f  ~rr  de  (4) 

Equation  (4)  taken  with  Equation  (1)  for  local  angle  of  attack,  when  computed  for  equi¬ 
valent  values  of  the  azimuth  angle,  gives  an  angle  of  attack  time  history  which  would  have  to  be 
imposed  on  the  wing  in  order  to  simulate  the  simultaneous  variation  of  both  aerodynamic  parameters. 

Figure  3  shows  wing  angle  of  attack  time  histories  determined  in  this  way  for  turbine  speed 
ratios  of  2,  3,  4  and  5.  For  the  speed  ratio  of  2  in  which  the  turbine  is  well  into  the  dynamic  stall 
regime,  and  the  maximum  angle  of  attack  is  30  degrees  and  the  distortion  of  the  cycle  is  considerable 
when  compared  with  a  simple  sine  wave.  Thus  as  the  blade  passes  from  the  advancing  to  retreating 
sides  of  the  rotor  azimuth,  a  at  first  increases  slowly  to  a  maximum  and  then  diminishes  rapidly 
through  zero  and  begins  the  negative  half  of  the  cycle.  When  maximum  negative  a  is  achieved  at  the 
same  rapid  rate  the  wing  returns  to  zero  at  the  initial  slow  rate. 

As  the  speed  ratio  increases,  maximum  a  decreases  as  also  does  the  large  difference  between 
the  slow  and  fast  portions  of  the  cycle. 

4.0  THE  OSCILLATING  AIRFOIL  EXPERIMENT 

The  experimental  program  consisted  of  the  aerodynamic  investigation  of  the  instantaneous 
pressure  and  dynamic  stall,  due  to  non-steady  motions  in  pitch,  of  a  NACA  0018  airfoil;  the  wing 
section  was  two-dimensional,  and  spanned  the  floor  and  the  ceiling  of  the  2m  X  3m  low  speed  wind 
tunnel.  The  angle  of  attack  time  history  was  periodic  and  approximated  the  azimuthal  aerodynamic 
flow  at  the  equatorial  plane  as  already  described  in  Section  3.0. 

Pressures  were  measured  on  both  the  model  and  wind  tunnel  walls,  and  the  time  history 
of  these  pressures  were  recorded  for  locations  along  the  chord  and  wall  belts.  The  wind  speed  which 
represented  the  equatorial  tip  speed  of  the  Magdalen  Islands  Turbine,  was  150  ft/sec  (45.7  m/sec). 
The  wing  frequency  was  constant  at  0.55  hz. 

4.1  Model 

The  wing  section  used  in  this  experiment  was  the  NACA  0018,  and  was  a  segment  of  the 
actual  aluminum  extruded  blade  profile  as  used  on  the  Magdalen  Islands  VAWT.  The  fidelity  of  the 
profile  to  the  standard  NACA  form  as  specified  in  Reference  (5),  was  not  exact.  The  quality  of  the 
surface  finish  was  also  below  the  usual  standard  for  wind  tunnel  models,  due  spanwise  waviness  and 
striations  resulting  from  the  extrusion  process. 

The  airfoil  was  mounted  in  a  vertical  plane  and  was  free  to  rotate  about  its  quarter-cord 
location.  The  lower  attachment  point  consisted  of  a  free  bearing,  located  beneath  the  tunnel  floor. 

The  upper  bearing  was  located  on  the  roof  turntable  and  supported  the  lever  mechanism  which  was 
connected  to  a  hydraulic  drive.  A  general  arrangement  of  the  wing  installation  is  shown  in  Figure  4. 

4.2  Hydraulic  Actuator 

The  hydraulic  actuator  and  control  system  was  designed  by  Davis  Engineering  Ltd.  of 
Ottawa,  and  was  used  as  a  driving  mechanism  for  the  oscillatory  airfoil  tests.  A  technical  description 
of  the  actuator  can  be  found  in  Reference  (6).  The  oscillatory  mechanism  consists  of  a  hydraulic 
servo  actuator,  an  electronic  controller  and  a  3000  psi,  2  gpm  hydraulic  power  pack. 


The  major  system  component  is  the  servo-actuator  which  oscillates  the  airfoil  through  a 
lever  arm  attached  to  a  shaft  on  the  wing.  This  device  is  mounted  on  a  bed-plate  which  also  carries 
the  associated  reservoirs,  valves  and  restrictors.  The  entire  assembly  is  fastened  to  the  upper  external 
wind  tunnel  structure.  The  linear  motion  of  the  actuator  is  converted  to  rotary  motion  by  means  of 
the  lever  arm.  The  position  of  the  wing  pinion  shaft  is  measured  by  a  rotary  transducer. 

The  servo-actuator  controls  the  position  of  the  airfoil  and  forces  the  hydraulic  actuator  to 
follow  an  externally-programmed  demand  signal.  This  is  accomplished  by  comparing  the  demanded 
wing  position  with  the  actual  position  as  sensed  by  the  rotary  transducer  mounted  in  the  shaft.  The 
resultant  error  signal  is  then  translated  into  a  drive  current  to  the  servo  valve;  this  in  turn  produces 
the  appropriate  hydraulic  flow  to  the  actuator,  thus  minimizing  the  wing  position  error. 

Under  the  dynamic  conditions  required  by  the  VAWT  simulation,  the  hydraulic  drive 
produced  an  airfoil  motion  which  was  within  1%  of  the  demanded  motion. 

4.3  Pressure  Measurements 

Surface  pressure  measurements  along  the  wing  chord  were  made  at  the  mid-span  position. 
These  locations  are  shown  schematically  in  Figure  5  and  in  Table  I.  Each  of  the  pressure  locations, 

45  in  all,  were  connected  to  the  header  of  a  Scanivalve  which  was  located  within  the  model  and 
rotated  with  it.  The  tubing  lengths  were  kept  to  a  length  of  about  2  to  3  feet  (61  to  91.5  cm). 

Since  the  dynamic  effects  of  the  tubing  on  fluctuating  pressures  is  considered  to  be  im¬ 
portant  in  measurements  such  as  this,  a  preliminary  investigation  was  made  as  to  possible  corrections 
to  the  pressures  which  might  be  necessary.  On  the  basis  of  Reference  (10)  and  also  as  a  result  of  using 
a  dynamic  signal  compensator  on  the  output  voltages,  it  was  concluded  that  there  were  no  serious 
errors  at  the  model  frequency  of  this  experiment. 

As  a  result  of  possible  unknown  dynamic  influences  due  to  the  wind  tunnel  walls,  pressures 
were  also  measured  on  each  side  wall  opposite  the  model.  The  location  of  the  wall  pressures  is  in¬ 
dicated  in  the  schematic  diagram  of  Figure  5,  and  also  in  Table  II.  The  wall  belt  pressures  were  also 
connected  to  a  Scanivalve. 

In  anticipation  of  significant  fluctuations  in  the  tunnel  flow,  both  tunnel  static  and  total 
pressures  were  measured  separately  using  a  standard  pitot  tube  placed  at  the  working  section  entrance 
plane.  Again,  due  to  expected  variations  in  tunnel  flow  all  pressures  were  referenced  to  atmospheric 
pressure  taken  outside  the  working  section.  This  pressure  Pr  is  on  the  reference  side  of  the  Scanivalve 
ansducers.  Thus  flow  dynamic  pressure,  referred  to  flow  static  pressure  is  determined  as  follows: 

PT  -P,  =  (Pt  -Pr)  ~  (P,-Pr)  (5) 

and  airfou  .  issure  coefficient  Cpa ,  also  defined  with  respect  to  flow  instantaneous  static  pressure  is: 

(Pa  -  Pr)  -  (Ps-Pr) 

Cpa  <PT  -  P,)  (6) 


Angle  of  attack  was  measured  relative  to  the  working  section  centre  line.  Positive  a  was 
defined  as  nose  right  (i.e.  toward  the  south  wall)  when  looking  upstream  from  above  (see  Fig.  4). 

4.4  Data  Acquisition  and  Reduction 

The  electrical  output  from  the  five  pressure  tranducers  and  the  rotary  transformer  was 
transcribed  on  to  six  channels  of  a  RACAL  FM  tape  recorder.  Since  the  Scanivalve  accesses  the 
pressures  sequentially,  it  was  necessary  to  record  for  each  of  the  wall  locations  (north  and  south)  and 
the  wing,  a  time  record  of  pressure  which  was  long  enough  to  encompass  several  oscillation  cycles. 
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(a)  WIND  TURBINE  POWER  CHARACTERISTICS 


(b)  TURBINE  INSTALLATION  -  MAGDALEN  ISLANDS 
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TABLE  II 

WALL  PRESSURE  LOCATIONS  (RELATIVE  TO  WING  QUARTER-CHORD) 

(NORTH  WALL) 

Port  No.  xw/c 


-2.58 

-1.92 

-1.59 

-1.26 

-0.93 

-0.60 

-0.27 

+0.06 

0.39 

0.72 

1.05 

1.38 

1.71 

2.04 

2.37 

2.70 

3.03 

3.36 

3.69 
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TABLE I 


I 


SCANIVALVE  PORT  NOS.  AND  WING  TAP  CO-ORDINATES 


Port  No. 


Wing  Co-Ordinate 
x/c  y/c 


Atmospheric  -  Reference 
Total  Pressure 


3 

0.946 

-0.0131 

4 

0.883 

-0.0226 

5 

0.821 

-0.0367 

6 

0.759 

-0.0462 

7 

0.698 

-0.0560 

8 

0.635 

-0.0642 

9 

0.573 

-0.0718 

10 

0.511 

-0.0780 

11 

0.463 

-0.0821 

12 

0.405 

-0.0860 

13 

0.343 

-0.0884 

14 

0.298 

-0.0900 

15 

0.256 

-0.0895 

16 

0.215 

-0.0868 

17 

0.174 

-0.0830 

18 

0.122 

-0.0740 

19 

0.098 

-0.0690 

20 

0.077 

-0.0635 

21 

0.057 

-0.0565 

22 

0.038 

-0.0480 

23 

0.0196 

-0.0350 

24 

0.0063 

-0.0170 

25 

0 

0 

26 

0.0063 

0.0170 

27 

0.0196 

0.0350 

28 

0.038 

0.0480 

29 

0.057 

0.0565 

30 

0.077 

0.0635 

31 

0.098 

0.0690 

32 

0.122 

0.0740 

33 

0.174 

0.0830 

34 

0.215 

0.0868 

35 

0.256 

0.0895 

36 

0.298 

0.0900 

37 

0.343 

0.0884 

38 

0.405 

0.0860 

39 

0.463 

0.0821 

40 

0.511 

0.0780 

41 

0.573 

0.0718 

42 

0.635 

0.0642 

43 

0.698 

0.0560 

44 

0.759 

0.0462 

45 

0.821 

0.0367 

46 

0.883 

0.0226 

47 

0"  46 

0.0131 
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3)  Integrated  values  of  the  wing  normal  and  chordwise  force  show  that  dynamic  stall  occurs 
much  later  than  for  steady  flow  resulting  in  an  increase  of  normal  force  of  about  20%  and 
an  increase  in  chordwise  force  of  about  40%,  when  the  wing  is  pitching  to  30  degrees. 

When  the  wing  was  pitched  to  maximum  angles  of  attack  well  beyond  the  steady  stalling 
angle,  a  considerable  degree  of  hysteresis  was  present.  This  characteristic  diminished  with 
maximum  angle  of  attack,  and  was  not  present  for  angles  of  attack  below  stall. 

4)  Observation  of  the  instantaneous  chordwise  and  thicknesswise  pressure  distributions  show 
how  the  growth  of  both  the  normal  and  chordwise  components  of  wing  force  originate, 
and  a  comparison  with  steady  flow  pressures  indicates  that  the  increase  of  leading  edge 
suction  pressure  is  about  50%  when  the  wing  is  pitching  to  a  maximum  angle  of  30  degrees. 
At  lower  maximum  pitch  angles,  particularly  below  steady  stall,  the  dynamic  effects  on 
suction  pressure  are  small. 

5)  The  data  were  not  corrected  for  blockage  or  flow  curvature  effects,  although  the  measure¬ 
ment  of  the  instantaneous  wall  pressures  were  obtained,  and  may  allow  further  investigation 
of  this  complex  subject. 

The  variation  of  flow  dynamic  pressure  when  the  wing  was  pitching  to  a  maximum  angle 
of  30  degrees,  was  about  6%  of  the  mean  value  which  itself  was  elevated  the  same  amount 
above  the  (so-called)  empty  tunnel  value.  This  elevation  and  variation  of  dynamic  pressure 
was  attributed  to  a  blockage  effect  probably  dynamic  in  orifin  and  which  diminished  with 
reduced  airfoil  maximum  pitch  angle. 
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At  smaller  maximum  angles  below  stall  the  pitching  moment  is  quite  small,  but  still  exhibits 
a  slight  anti-symmetry  with  larger  values  in  the  negative  angle  range.  The  large  negative  moments  which 
occur  after  stall  have  been  attributed  to  free  vortices  which  have  been  shed  regularly  from  the  leading 
edge  and  pass  over  the  aft  portions  of  the  chord.  This  increases  local  loading  and  causes  the  wing  to 
pitch  down.  The  vortices  are  also  related  to  the  augmentation  of  both  lift  and  drag  during  the  dynamic 
stall  phase.  These  large  negative  moments  are  also  present  under  steady  flow  conditions,  as  shown  in 
Figure  25a. 

6.6  Comparison  between  Steady  and  Unsteady  Airfoil  Characteristics 

Measurements  of  airfoil  section  pressures  under  steady  flow  conditions  were  made  in  a 
separate  experiment.  Figure  26a  shows  at  an  angle  of  attack  of  20  degrees,  a  comparison  between  a 
steady  chordwise  pressure  distribution  and  a  dynamic  chordwise  pressure  distribution  when  the  airfoil 
was  pitching  to  amax  =  30  degrees.  The  increase  of  maximum  suction  pressure  at  the  nose  was  from 
4.0  to  7.0.  Similarly  the  positive  loop  in  Figure  26b  indicates  a  considerable  increase  in  the  upwind 
suction  force  at  this  incidence. 

Similar  graphs  for  the  unstalled  airfoil  (amax  =  12)  show  a  much  more  modest  increase  in 
suction  pressure  at  due  to  dynamic  effects  at  the  incidence  for  maximum  lift  (Fig.  27). 

The  integrated  values  of  Cn  and  Ctmay  also  be  compared.  The  steady  value  of  Cn  is  not  too 
different  from  the  corresponding  dynamic  value  except  after  stall  occurs.  The  dynamic  effects  seem 
to  produce  a  much  larger  hysteresis  and  asymmetry.  In  the  case  of  the  leading  edge  force  Ct  the  dif¬ 
ference  between  static  and  dynamic  load  is  larger  and  the  asymmetry  between  positive  and  negative 
angular  ranges  is  much  more  pronounced  (Fig.  28). 

7.0  CONCLUSIONS 

The  following  conclusions  can  be  made  regarding  the  data  presented  for  the  NACA  0018 
airfoil  oscillating  in  pitch  to  various  maximum  angles  of  attack.  The  purpose  of  the  test  was  two-fold: 
(a)  to  simulate  the  dynamic  aerodynamic  behaviour  of  a  VAWT  blade  section  at  the  equator  in  a 
quasi-two-dimensional  flow  environment,  and  (b)  to  investigate  the  characteristics  of  the  airfoil  under 
dynamic  stall  conditions  which  occur  at  turbine  speed  ratios  (Vt/V)  of  3  or  less. 

Measurements  were  made  of  wing  surface  pressures  and  flow  static  and  total  pressures, 
relative  to  a  fixed  atmospheric  reference.  The  time  histories  of  one  complete  cycle  of  oscillation  were 
obtained  from  the  average  of  five  cycles  of  data.  Pressure  coefficients  were  computed  and  integrated 
to  give  normal  and  chordwise  force  coefficients. 

1)  The  pressure  time  histories  on  the  wing  were  characterized  by  intervals  of  high  suction 
pressure  which  occurred  twice  per  cycle  near  the  nose  as  the  wing  oscillated  from  side  to 
side.  These  twin  pressure  peaks  which  occurred  on  both  sides  of  the  airfoil  diminished  to  a 
single  configuration  as  the  chordwise  stations  further  rearward  experience  suction  only  once 
per  cycle. 

2)  Dynamic  stall  occurred  near  an  angle  of  attack  of  20  degrees  when  the  wing  was  pitching 
to  amax  =  30  degrees.  The  normal  stalling  angle  under  steady  flow  conditions  at  the  same 
Reynolds  number  was  15  degrees.  Nose  suction  hung  on  for  some  time  after  flow  breakdown 
and  complete  separation  did  not  occur  even  at  30  degrees  incidence  as  evidenced  by  the 
high  suction  peaks  at  the  nose.  There  was  no  clear  evidence  from  the  wing  pressure  time 
history  carpets,  of  a  pressure  wave  induced  by  a  free  vortex  travelling  downwind  over  the 
chord  of  the  wing. 


At  angles  of  attack,  which  under  steady  flow  conditions  would  not  produce  stalling,  the 
wing  pressures  are  not  significantly  increased  by  the  dynamic  motion  of  the  airfoil  although 
the  suction  peaks  at  the  nose  were  maintained. 


At  zero  incidence  the  area  of  the  positive  loops,  which  produce  an  upwind  force,  approxi¬ 
mately  balance  the  negative  loop  of  positive  pressures  which  act  at  the  nose,  and  result  in  a  downwind 
or  drag  force.  The  net  drag  at  zero  incidence  is  determined  by  this  balance  of  pressures  and  the  skin 
friction  drag  which  is  normally  determined  by  wake  traverse. 

As  incidence  increases  the  area  of  the  positive  pressure  loop,  composed  of  those  pressures 
originating  on  the  suction  side  of  the  wing,  increases  steadily,  and  result  in  a  large  upwind  pressure 
force.  This  force  is  still  active  after  the  wing  has  stalled  due  to  the  fact  that  the  oncoming  flow  which 
has  impinged  on  the  leading  edge  and  is  still  able  to  negotiate  the  nose  curvature  without  separating. 
This  phenomenon  is  considerably  enhanced  under  dynamic  conditions  and  the  behaviour  of  the 
unsteady  leading  edge  boundary  layer  is  presumably  a  key  factor  in  the  understanding  of  airfoil 
dynamic  stall.  At  a  maximum  incidence  of  30  degrees,  the  suction  pressures  are  still  active  and  a  net 
upwind  force  exists,  although  much  diminished. 

When  the  airfoil  pitches  to  a  smaller  maximum  incidence,  corresponding  to  higher  speed 
ratios,  it  does  not  experience  a  complete  stall  and  the  dynamic  characteristics  were  not  noticeably 
different  from  those  of  steady  flow.  The  leading  edge  pressures  embodied  in  the  positive  loop  of  the 
pressure  diagram  produce  an  upwind  force  which  also  prevails  after  flow  breakdown  has  occurred  on 
the  airfoil  (Figs.  15  to  18). 

6.5  Airfoil  Normal  and  Chord  Force  Coefficients 

Figures  19  to  22  show  the  behaviour  of  airfoil  normal  and  chord  force  coefficients  derived 
from  integration  of  the  instantaneous  chordwise  and  thicknesswise  pressures.  The  maximum  values 
of  Cn  and  Ct  depend  upon  the  maximum  angle  to  which  the  wing  pitches.  For  a  speed  ratio  of 
2(amax  =  30),  Cnmax  occurs,  nears  18  degrees  and  has  a  value  of  1.38.  The  corresponding  value  of 
Ctmax  is  0-4.  Both  values  of  Cn  and  Ct  are  considerably  higher  than  the  equivalent  steady  values,  due 
to  dynamic  stall  effects. 

The  normal  and  chord  force  coefficients  shown  in  Figure  19  exhibit  characteristics  of  both 
hystereses  and  non-symmetry  which  are  not  present  in  the  same  form  in  the  steady  airfoil  character¬ 
istics.  There  is  no  satisfactory  explanation  for  this  vehaviour  at  this  time.  The  normal  and  chord  force 
characteristics  for  smaller  maximum  angles  of  pitch  (higher  speed  ratios)  produce  lower  values  of  Cn 
and  C, ,  and  do  not  exhibit  the  significant  hysteresis  and  non-symmetry  of  the  stalled  airfoil. 

Figures  23  and  24  show  the  airfoil  instantaneous  lift  and  drag  coefficients  derived  from  the 
integrated  pressure  forces  and  plotted  against  incidence.  The  wing  achieved  a  maximum  Cj  of  1.4  at 
an  angle  of  attack  of  20  degrees  (amax  =  30).  This  may  be  compared  favourably  with  a  steady  value 
of  1.16  in  which  the  wing  stalls  at  15  degrees.  At  a  lower  value  of  maximum  incidence  (amax  =  12, 
Fig.  23b),  Clmax  =  1,  and  the  lift  curve  exhibits  almost  no  hysteresis.  The  lift  slope  for  a  maximum 
angle  of  12  degrees  was  0.08  per  degree.  The  lift  slope  for  amax  =  30  in  which  dynamic  stall  occurred, 
was  on  the  average,  less  than  for  the  unstalled  case  pitching  to  a  maximum  angle  of  12  degrees. 

The  drag  coefficient  Cd  is  plotted  in  Figure  24  for  maximum  pitch  angles  above  and  below 
stall  (amax  =  30  and  12).  Drag  increases  rapidly  with  incidence  and  when  the  airfoil  is  stalled,  but  has 
a  lower  value  on  the  fast  return  portion  of  the  pitch  cycle,  thus  corresponding  to  the  hysteresis  in  lift. 
Drag  at  zero  incidence  is  not  believed  to  be  accurately  determined  in  this  experiment,  however  a  value 
as  determined  from  wake  measurements  with  the  wing  fixed  is  Cdo  =  0.01. 

Wing  pitching  moment  coefficient  (about  wing  quarter-chord),  is  plotted  in  Figure  25  for 
maximum  wing  pitch  angles  of  30  and  12  degrees  respectively.  Positive  moment  is  in  the  direction  of 
increasing  angle  of  attack.  The  pitching  moment  curve  in  Figure  25a,  a  =  30  degrees,  exhibits  neutral 
stability  between  plus  and  minus  15  degrees.  When  the  wing  is  stalled  and  at  angles  greater  than 
25  degrees,  a  large  negative  moment  exists  and  some  hysteresis  is  present.  The  general  shape  of  the 
moment  characteristic  is  anti-symmetric,  but  less  so  in  the  negative  angle  range  (see  also  the  Cn  and 
Ct  curves,  Fig.  19). 


Figure  7  shows  pressure  distributions  for  Vt/V  =  2  (amax  =  30).  Under  static  conditions 
with  the  wing  fixed  at  30  degrees,  the  flow  would  be  completely  separated,  with  high  negative  pres¬ 
sures  appearing  only  over  a  small  portion  of  the  rounded  nose.  This  figures  shows,  however  that  the 
negative  peak  pressures  are  much  larger  in  magnitude  under  dynamic  conditions  than  those  for  fixed 
incidence,  and  that  these  pressures  are  maintained  at  a  high  level  well  aft  of  the  nose  until  close  to 
maximum  angle  of  attack.  This  suggests  that  flow  breakdown  and  stall  is  delayed  under  dynamic 
conditions.  Surface  pressures  on  the  rear  portions  of  the  airfoil  after  flow  breakdown  has  occurred 
correspond  to  ambient  static  pressure  although  their  time  histories  indicate  a  negative  pressure  peak 
close  to  maximum  angle  of  attack.  Without  an  extensive  flow  investigation  it  is  difficult  to  explain 
this  suction  pressure  on  the  rear  of  the  airfoil;  the  effect  of  separated  flow  on  the  end  walls,  could  be 
important  in  this  regard,  producing  strong  vortex-induced  velocities  acting  on  the  rear  of  the  airfoil. 

The  time  history  of  pressures  near  the  nose  of  the  airfoil  during  the  slow  portion  of  the 
cycle,  indicates  a  rapid  rise  in  negative  pressure  which  peaks  just  before  amax.  Beyond  amax  pressure 
falls  rapidly,  as  angle  of  attack  decreases  during  the  fast  portion  of  the  cycle  and  flow  recovery  is 
re-established  on  the  suction  side  of  the  airfoil.  As  the  wing  incidence  passes  through  zero  and  be¬ 
comes  negative,  the  pressures  reverse  sign  reaching  a  maximum  positive  value  of  1.0  as  stagnation 
pressure  is  encountered. 

It  can  be  seen  from  the  time  histories  that  pressures  near  the  nose  experience  two  maxima 
in  both  suction  and  positive  pressure  levels  as  they  pass  alternately  through  regions  of  maximum 
surface  velocity  and  stagnation  respectively. 

As  the  maximum  angle  of  attack  decreases  (increasing  speed  ratio),  the  pressure  peaks  near 
the  nose  of  the  airfoil  gradually  decrease  in  magnitude  but  still  exhibit  the  double  maxima  as  the 
stagnation  point  oscillates  from  side  to  side.  Suction  pressure  is  also  maintained  further  along  the 
wing  chord,  indicating  that  the  flow  separation  is  confined  mainly  to  the  trailing  edge  region.  In 
Figure  10  for  example,  amax  =  12,  the  variation  of  suction  pressure  on  the  forward  portion  of  the 
airfoil  reaches  maximum  values  of  4.0;  a  value  close  to  that  of  the  fixed  airfoil. 

6.4  Airfoil  Instantaneous  Surface  Pressure  Distributions 

Figures  11  and  12  show  the  chordwise  and  thicknesswise  pressure  distributions  on  the 
surface  of  the  airfoil  at  a  number  of  angle  of  attack  positions  during  a  complete  cycle  of  oscillation. 

In  Figure  11,  V,/V  =  2,  amax  =  30,  the  pressure  distributions  show  the  development  of  leading  edge 
suction,  lift  and  flow  breakdown  as  the  airfoil  pitches  slowly  to  maximum  positive  incidence,  and 
then  returns  quickly  through  zero  to  maximum  negative  incidence,  followed  by  slow  return  to  zero. 

The  chordwise  pressures  represent  aerodynamic  flows  at  all  of  the  speed  ratios  and  illustrate 
the  effect  of  maximum  pitch  incidence  on  growth  of  lift  and,  in  particular,  the  development  of  flow 
separation  and  reattachment. 

Pressure  coefficients  plotted  across  the  thickness  of  the  airfoil  are  shown  in  Figure  12. 
These  graphs  illustrate  clearly  the  origin  of  the  leading  suction  force  and  its  growth  with  incidence. 
The  magnitude  of  the  net  force  acting  along  the  wing  chord  is  determined  from  the  differences  be¬ 
tween  the  positive  and  negative  loops  in  the  pressure  diagrams  of  Figure  12.  This  is  illustrated  dia- 
grammatically  as  follows. 


swassss  msea»r  sassm  msem  mm 


-7- 


6.0  EXPERIMENTAL  RESULTS 

6.1  Wind  Tunnel  Test  Program 

The  wind  tunnel  test  program  of  the  Oscillating  Airfoil  was  run  at  a  constant  nominal 
dynamic  pressure  of  27  psf;  and  a  constant  pitch  frequency  of  0.55  hz.  A  slow  pitch  rate  was  also 
done  in  order  to  furnish  steady  aerodynamic  characteristics  of  the  airfoil. 

The  simulation  of  the  dynamic  variation  of  angle  of  attack  and  resultant  velocity  time 
histories  was  done  for  turbine  speed  ratios  of  2,  3,  4  and  5.  These  values  correspond  to  maximum 
angles  of  attack  above  and  below  stall,  and  represent  the  various  performance  regions  of  VAWT 
operation  (Fig.  1). 

Measured  data  is  presented  in  the  form  of  pressure  coefficients  Cp  at  various  locations  on 
the  wing  and  walls.  These  coefficients  show  the  time  history  of  the  unsteady  pressures  as  the  wing 
pitches  through  its  angle  of  attack  program.  Also  constructed  are  the  conventional  chordwise  and 
thicknesswise  pressure  distributions  throughout  the  pitch  cycle.  These  pressures  when  integrated, 
furnish  normal  and  chord  force  coefficients  under  dynamic  aerodynamic  flow  conditions. 

The  data  as  presented  in  this  report  is  uncorrected  for  any  of  the  usual  blockage  and  flow 
curvature  effects.  The  discussion  of  the  application  of  conventional  wind  tunnel  corrections  in  un¬ 
steady  type  of  airfoil  experiments  similar  to  the  present  case  can  be  found  in  Reference  (7).  It  was 
suggested  that  although  this  type  of  corrective  method  can  be  satisfactorily  applied  to  lift,  it  produces 
large  errors  in  drag  due  to  the  difficulty  in  resolving  the  difference  between  two  large  quantities; 
namely  Cn  and  Ct  particularly  near  maximum  lift. 

In  this  experiment  the  wall  pressures  were  measured  at  the  same  time  as  the  wing  pressures, 
so  that  the  instantaneous  influence  of  the  walls  could,  in  principle,  be  determined  from  wall  signature 
techniques  for  steady  flow  such  as  those  described  in  References  (8)  and  (9).  Unsteady  wall  effects 
even  at  the  low  frequencies  of  this  experiment  may  not  be  quasi-steady  however,  and  secondary  flows 
may  require  further  correction.  Accordingly,  it  was  thought  that  a  complete  analysis  of  the  problem 
would  be  outside  the  scope  of  this  report  but  could  be  addressed  in  a  future  investigation. 

6.2  Time  Variation  of  Tunnel  Static  and  Dynamic  Pressure 

The  time  variation  of  tunnel  dynamic  and  static  pressures  in  the  flow  at  the  working  section 
entrance  plane  is  presented  in  Figure  6  for  one  cycle  of  oscillation.  For  a  wing  oscillation  of  30  deg¬ 
rees  corresponding  to  a  turbine  speed  ratio  of  2,  and  a  wing  angle  of  attack  well  beyond  stall  (a  max  = 
30)  the  peak  to  peak  variation  of  dynamic  pressure  was  6%  of  the  nominal  reference  value  of  27  psf. 
The  mean  value  of  dynamic  pressure  appears  to  be  elevated  about  6%  above  the  reference  value.  The 
corresponding  variation  of  static  pressure  was  6%  of  tunnel  dynamic  pressure,  and  its  mean  value  was 
also  elevated  well  above  the  reference  value.  For  a  maximum  angle  of  attack  of  12  degrees  (Vt/V  =  5), 
variations  in  both  tunnel  dynamic  and  static  pressures  are  small  and  their  mean  values  are  close  to  the 
reference  pressures.  These  pressure  variations  and  their  fluctuations,  particularly  at  angles  of  attack 
beyond  stall,  are  probably  a  measure  of  blockage  resulting  from  both  conventional  sources  and  oscil¬ 
lation  effects  on  the  tunnel  average  flow. 

6.3  Pressure  Time  Histories  on  the  Wing 

The  time  histories  of  surface  pressures  on  the  wing,  Cpo ,  are  presented  in  the  form  of  carpet 
plots  for  the  four  turbine  speed  ratios  (Vt/V  =  2,  3,  4,  5)  which  were  simulated  in  this  experiment. 

In  this  form  of  presentation,  (Figs.  7  to  10),  both  surfaces  of  the  symmetrical  airfoil  pressures  (i.e. 
north  and  south)  are  presented  and  the  time  histories  for  each  port  are  offset  so  as  to  illustrate  the 
progression  of  unsteady  aerodynamic  effects  along  the  airfoil  chord.  The  time  base  represents  one 
complete  cycle  of  oscillation. 
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The  angle  of  attack  channel  6,  was  used  to  correlate  all  other  channels  with  a  known  time 
base,  since  its  amplitude  and  frequency  were  controlled  and  unchanging. 

5.0  CALCULATION  OP  PRESSURE  COEFFICIENTS 
Wall  pressure  (south) 

CpWS(t)  =  (Kl*(Ej  -  E100)  -  (ps  -  pr))/(pt  -  p,) 

Wall  pressure  (north) 

Cpwnft)  =  (K2*(E2  -  E200)  -  (ps  -  pr))/(pt  -  ps) 

Airfoil  pressure 

Cpa(t)  =  (K3*(Ej  -  E300)  -  (pg  -  pr))/(pt  -  ps) 
where  the  static  pressure  is 

(Ps  “  Pr)  =  K5*(ES  -  E500) 
and  total  pressure  is 

(Pt  ~  Pr)  =  K6*(E6  -  E600) 
and  the  instantaneous  flow  dynamic  pressure  is 

(Pt  -  P,)  *  (Pt  -  Pr)  ~  (Ps  -  Pr) 

5 .1  Integration  of  Wing  Pressure  Coefficients 

The  wing  pressure  coefficients  for  each  of  the  45  ports  of  the  Scanivalve  are  integrated  over 
the  chordwise  and  thickness  dimensions  to  give  coefficients  of  normal  force  C„ ,  chord  force  Ct  and 
pitching  moment  Cm .  (Wing  co-ordinates  are  listed  in  Table  I). 

Normal  Force 

C"  =  i=3  Cpa*dx"  ?6  CP#*dX 
Chord  Force 

Ct  =  i?3  CP»*dy"  J  cpa*dy+  2  Cpa*dy 
Pitching  Moment  (1/4  chord) 

Cmc/4  *  2  (0.25  -  x)*Cpg*dx  -  2  (0.25  -  x)*Cpa*dx 
i“  3  26 

5.2  Wing  Instantaneous  Lift  and  Drag  Coefficients 

Cj  *  C„  *Cosa  +  Ct  *Sina 


Cd  ■  -(Cn*Sina  -  C,*Cosa) 


Tunnel  static  pressure,  total  pressure  and  wing  angle  were  recorded  simultaneously  and  continuously. 
Thus  for  each  of  45  wing  pressure  locations,  and  21  wall  pressure  locations,  a  time  record  was 
produced  which  could  be  correlated  with  both  tunnel  flow  conditions  and  wing  angle. 


It  was  necessary  to  calibrate  each  pressure  channel,  and  this  was  done  by  connecting  the 
working  section  total  pressure  as  a  reference  value  to  port  #01  on  each  Scanivalve.  With  the  wing 
motionless  and  set  at  zero  incidence,  the  tunnel  was  run  at  a  nominal  dynamic  pressure  of  27  psf 
(corresponding  to  150  ft/sec  (45.7  m/sec)),  and  the  pressure  signal  was  recorded  on  all  five  pressure 
channels.  The  actual  total  and  static  pressures  relative  to  atmosphere,  were  measured  independently 
using  a  Betz  Manometer.  Their  values,  for  a  tunnel  pressure  of  27  psf  were  28.146  psf  and  1.091  psf 
respectively. 

Angle  of  attack  was  calibrated  in  a  similar  way  by  setting  the  wing  to  a  known  value  and 
recording  the  signal  transmitted  by  the  rotary  transducer  on  channel  6. 

Following  the  conversion  of  the  RACAL  analogue  signals  to  digital  form,  the  data  were 
processed  to  furnish  the  desired  pressure  and  force  coefficients. 

4.5  Calculation  of  Pressure  Sensitivities 

Calibration  of  each  channel  transducer  took  place  at  the  beginning  of  each  run  with  the 
wing  set  at  zero  degrees.  The  differences  between  wind  on  and  wind  off  portions  of  the  time  record 
were  related  to  the  separately  measured  tunnel  total  and  static  pressures  already  specified.  Calculation 
of  pressure  sensitivities  for  each  channel  was  as  follows: 

Channel  1  —  Wall  pressure  (south) 

K1  =  28.146/(E101  -E100) 

Channel  2  —  Wall  pressure  (north) 

K2  =  28.146/(E2oi  -  E200) 

Channel  3  —  Airfoil  pressure 

K3  =  28.146/(E301  -  E300) 

Channel  4  —  Total  pressure 

K4  =  28.146/(E401  -  E400) 

Channel  5  —  Static  pressure 

K5  =  1.091/(E501  -Esoo) 

where  Kl,  K2,  K3,  K4,  K5  are  the  respective  pressure  sensitivities.  The  corresponding  values  of  En0, 
and  En00  are  the  wind  off  and  wind  on  output  voltage  magnitudes,  with  the  wing  set  to  zero  incidence. 

The  angle  of  attack  sensitivity  K6  was  determined  in  a  similar  way. 

4.6  Sampling  and  Averaging  of  the  Data 

During  the  analogue  to  digital  processing  the  sampling  rate  was  250  per  second  and  about 
20  seconds  of  information,  corresponding  to  10  wing  oscillations,  was  taken  for  each  Scanivalve  port 
position.  After  editing,  five  complete  cycles  in  the  middle  of  the  time  record  were  averaged  to  mini¬ 
mize  any  random  fluctuations  which  might  be  present. 
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STATIC  PRESSURE  LOCATIONS  ON  THE  WING  AND  WALLS 
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FIG.  6:  WIND  TUNNEL  DYNAMIC  AND  STATIC  PRESSURE  TIME  HISTORIES 

DURING  ONE  CYCLE 


FIG.  7:  WING  PRESSURE  TIME  HISTORIES,  a_  =  ±30°,  —  =  2,  f  =  0.55  hz 


FIG.  8:  WING  PRESSURE  TIME  HISTORIES,  a_  =  ±20°,  ^  =  3,  f  =  0.55  hz 
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FIG.  11b:  WING  CHORDWISE  PRESSURES  DURING  NEGATIVE  HALF-CYCLE, 

On,  =  ±30°,  ^  =  2,  f  -  0.55  hz 
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FIG  12a:  WING  THICKNESSWISE  PRESSURES  DURING  POSITIVE  HALF-CYCLE, 

DO 

o„,  =  ±30°,  ^  =  2,  f  =  0.55  hz  (Cont’d) 
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FIG.  14a:  WING  THICKNESSWISE  PRESSURES  DURING  POSITIVE  HALF  CYCLE, 

DO 

o_  =  ±20°.  ^  -  3.  f  -  0.55  hz  (Cont'd) 
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FIG.  12b:  WING  THICKNESSWISE  PRESSURES  DURING  NEGATIVE  HALF-CYCLE, 
a„,  =  ±30°,  —  =  2,  f  =  0.55  hz  (Cont'd) 
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FIG.  18b:  WING  THICKNESSWISE  PRESSURES  DURING  NEGATIVE  HALF-CYCLE, 
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(b)  CHORD  FORCE  Cy 


FIG.  19:  AIRFOIL  FORCE  COEFFICIENTS  CN,  CT,  DERIVED  FROM  PRESSURE 
DISTRIBUTIONS  =  2.0,amax  =  ±30°,  f  =  0.55  hz 
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FIG.  21:  AIRFOIL  FORCE  COEFFICIENTS  CN,  Cj,  DERIVEO  FROM  PRESSURE 
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SUMMARY/SO  MMAIRE 

">  This  report  presents  aerodynamic  data  on  an  NACA  0018  airfoil  oscillating  in  pitch  in  the/\ 
2m  X  3m  Low  Speed  Tunnel.  The  purpose  of  the  test  was  to  simulate  the  dynamic  aerodynamic 
behaviour  of  a  Vertical  Axis  Wind  Turbine  blade  section  at  the  equatorial  plane  under  the  dynamic 
stall  conditions  which  may  occur  at  low  ratios  of  tip  speed  to  wind  speed. 

Measurements  were  made  of  wing  surface  pressure  time  histories  at  various  chordwise 
locations,  during  a  complete  cycle  of  oscillation.  The  dynamic  effects  on  nose  suction  pressures  were 
seen  to  increase  their  magnitude  and  to  delay  flow  breakdown  to  a  higher  angle  of  attack. 

Integration  of  wing  pressures  furnished  normal  and  chordwise  forces  during  the  angle  of 
attack  cycle.  The  behaviour  of  these  parameters  showed  that  dynamic  stall  for  wing  angles  of  attack 
up  to  30  degrees,  occurred  about  5  degrees  later  than  for  the  equivalent  steady  flow  case.  This  phe¬ 
nomenon  resulted  in  an  increase  in  normal  force  of  about  20%  and  an  increase  in  chordwise  force  of 
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about  40%  when  the  wing  is  pitching  to  30  degrees. 
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